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ABSTRACT 

We present a high resolution, SparsePak two-dimensional velocity field for the center of the low surface bright- 
ness (LSB) galaxy DDO 39. These data are a significant improvement on previous HI or Ha long slit data, yet the 
inner rotation curve is still uncertain due to significant noncircular and random motions. These intrinsic uncertain- 
ties, probably present in other LSB galaxies too, result in a wide range of inner slopes being consistent with the 
data, including those expected in cold dark matter (CDM) simulations. The halo concentration parameter provides 
a more useful test of cosmological models than the inner slope as it is more tightly constrained by observations. 
DDO 39's concentration parameter is consistent with, but on the low end of the distribution predicted by CDM. 

Subject headings: galaxies: dwarfs — galaxies: halos — galaxies: kinematics and dynamics 



velocity field of the LSB galaxy DDO 39. This galaxy was 
selected because of its suitable inclination of 49°, low surface 
brightness = 24.4 mag arcsec" 2 , S waters & Balcells 2002), 
disk dominated radial surface brightness profile with a disk 
scale length of h = 3.5 kpc at an adopted distance of 12.8 Mpc 
(S99), regular optical and HI morphology, and well-behaved 
HI kinematics (see Fig. 1). These properties make DDO 39 an 
ideal galaxy to study the kinematics in the central regions of an 
LSB galaxy. 

2. OBSERVATIONS AND DATA REDUCTION 

DDO 39 was observed with the SparsePak integral field 
unit on the 3.5m WIYN telescope on January 11, 12 and 13, 
2002. SparsePak is a fiber-optic array containing 82 fibers, each 
4.7" in diameter and separated by 5.6", arranged in a sparsely 
packed grid, 72" on each side, with a small, nearly-integral 
core. Seven sky fibers are spaced between 60" and 90" away 
from the central fiber. The grid can be fully sampled in three 
pointings. For a detailed description, see Bershady et al. (in 
preparation). The fibers are fed to the WIYN Bench Spectro- 
graph. We used the 860 1/mm grating in second order at 6600 
A, yielding a FWHM velocity resolution of 66 km s" 1 . In the 
first two nights we filled the grid with three 3 x 20 min integra- 
tions; in the third night we integrated 3x1 hour on the central 
pointing. 

The IRAF dohydra package was used for data reduction and 
extraction of the spectra. Sky subtraction was done by creat- 
ing a two-dimensional image of extracted spectra, filtering out 
the Ha lines, and subtracting a low order baseline (Bershady et 
al, in preparation). To determine the actual telescope pointing 
and offsets, we calculated the expected relative continuum flux 
for all fibers from an /?-band image of DDO 39, and compared 
this to our observations. By minimizing the x 2 difference, the 
pointing offsets could be determined with an accuracy of ~ 0.5" 
(Swaters et al., in preparation). 

We constructed a sparse velocity field by fitting Gaussians to 
spectra with Ha emission stronger than three times the noise 
level and placing the derived velocities in a map at the corre- 
sponding positions. To construct a contiguous velocity field 
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1. INTRODUCTION 

Because low surface brightness (LSB) galaxies are most 
likely dominated by dark matter at all radii (e.g., de Blok & Mc- 
Gaugh 1997; Verheijen 1997; Swaters, Madore, & Trewhella 
2000), they are ideal for studying the power law slopes of the 
central dark matter density distributions p{r) oc p~ a . Cosmo- 
logical simulations indicate a depends on the nature of the dark 
matter (e.g., Navarro, Frenk, & White 1997, hereafter NFW; 
Fukushige & Makino 1997; Moore et al. 1999; Colin, Avila- 
Reese, & Valenzuela 2000; Dave et al. 2001; Knebe et al. 2001). 
Thus, a measurement of a may provide constraints on the na- 
ture of dark matter in galaxies and theories of galaxy formation. 

Unfortunately, it is difficult to measure the inner slope a ob- 
servationally. HI observations have relatively low angular reso- 
lution and are affected by beam smearing (Swaters 1999, here- 
after S99; Swaters et al. 2000; van den Bosch et al. 2000; Mc- 
Gaugh, Rubin, & de Blok 2001); hence HI observations gen- 
erally are less suited to measure a accurately. Ha long-slit ob- 
servations provide higher angular resolution, but studies based 
on such observations find conflicting results. Some find LSB 
galaxies are consistent with steep inner slopes (Swaters 2001; 
Swaters et al. 2003, hereafter SMvdBB), yet others find them 
inconsistent with steep slopes (de Blok et al. 2001a,b; de Blok 
& Bosma 2002, hereafter dBB; Marchesini et al. 2002). This 
inconsistency remains even when identical datasets are used 
(SMvdBB). The apparent discrepancy can be explained, in part, 
by systematic effects of slit width, seeing, slit offsets and galaxy 
inclination, which all lead to an underestimate of a (SMvdBB). 

To avoid systematic effects of long-slit observations, and 
to map possible non-circular motions, it is necessary to ob- 
tain high spatial resolution, two-dimensional velocity fields 
(see also Beauvais & Bothun 1999). Blais-Ouellette et al. 
(2001) presented such observations for the dwarf LSB galaxies 
NGC 3109 and IC 2574. Unfortunately, both galaxies are less 
suited to address the core issue; NGC 3109 is close to edge-on, 
and IC 2574 has a poorly sampled velocity field and a perturbed 
interstellar medium (Walter & Brinks 1999). 

In this Letter, we present a high-resolution, two-dimensional 
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FIG. 1. — Hi velocity field from Swaters et al. (2002) (left panel), SparsePak 10" velocity field (middle panel), and SparsePak 5.5" velocity field (right panel) 
superposed on the optical image from Swaters & Balcells (2002). Steps between isovelocity contours are 10 km s _1 in the left two panels, and 5 km s _1 in the right 
panel. The dotted lines represent the approaching side, the first full line indicates the systemic velocity. The beam size is given in the lower right of each panel, the 
cross indicates the galaxy center. Most bright objects are foreground stars. 



for graphical presentation, the points in the sparse velocity field 
were interpolated after weighting with a Gaussian beam of 5.5" 
and 10". These velocity fields, shown in Fig. 1, highlight the 
kinematics in the central regions and over the entire SparsePak 
field of view. Because of the interpolation, the isovelocity con- 
tours may be uncertain, especially at the edges. 

3. RESULTS 
3.1. The rotation curve 

To determine the galaxy orientation parameters, we fitted a 
tilted ring model to the sparse velocity field. Because of the ap- 
proximately linear rise of the rotation curve (RC), it was not 
possible to determine the center or the inclination kinemati- 
cally from the SparsePak data. Instead, we used the optical 
center at 5 ft 56 d 37 .5 s 75°18'56" (J2000). Although the inner 
rotation curve shape depends somewhat on the choice of cen- 
ter, the main results of this Letter do not. For the inclination we 
used 49°, derived from the HI observations. The position angle 
and systemic velocity were found to be consistent with the 133° 
and 818 km s" 1 of the HI observations, and we used the latter. 

With these orientation parameters, all points in the sparse 
velocity field were corrected for projection effects. To deter- 
mine the rotation velocity, points were averaged in concentric 
annuli with a width of 5". Each point was assigned a weight 
of (cos</>/cr ra d) 2 , where <j) is the angle in the plane of the galaxy 
with respect to the major axis, and cr ra d is the error on the radial 
velocity, with an imposed minimum of 4 km s" 1 . The error on 
the rotation velocity was taken to be the quadratic sum of the 
formal error and one fourth of the difference between the ap- 
proaching and receding sides, with an imposed minimum error 
of 2 km s -1 . The RC derived in this way was combined with the 
HI RC from S99 to create the hybrid RC shown in Fig. 2a. 

Fig. 2a shows good agreement between our RC and the RC 
derived by S99 from HI data. Fig. 2b compares dBB's long-slit 
data and our data within 2.5" of the major axis. There is good 
general agreement, although the velocities derived by dBB tend 
to be somewhat closer to the systemic velocity. Fig. 2c com- 
pares our RC to the one derived by dBB. Although consistent 
within the errors, their RC is substantially lower where inde- 
pendent optical data are used (dBB also used S99's data beyond 
~ 90"). The reason for this is unclear, but the fact that we used 
two-dimensional data and dBB spline-interpolated their binned 



data to derive their RC are both likely to play a role. 



3.2. Noncircular motions 

Above we have assumed that the gas moves on circular or- 
bits. However, the Ha velocity fields in Fig. 1 reveal signifi- 
cant noncircular motions. Because the velocity fields might be 
affected by interpolation, we have also plotted the data directly. 
In Fig. 3 we plot the observed radial velocities with respect to 
systemic velocity in different radial intervals, normalized by the 
rotation velocity. At large radii the radial velocities follow the 
curve expected for circular rotation closely. On the other hand, 
in the innermost regions there are clear deviations from simple 
circular motions. 

Fig. 3 also shows that there is significant scatter around the 
best fit rotation velocity. The scatter (before normalization) is 
around 9 km s" 1 , independent of radius. The scatter is in part 
due to observational uncertainties. Measurement errors depend 
on the signal-to-noise ratio and range from 2 to around 8 km 
s" 1 for the weakest lines, and average about 4 km s" 1 . Addi- 
tional scatter may arise if the intensity weighted position of the 
Ha over the fiber "beam" does not coincide with the center of 
that fiber. The maximum contribution of this effect, modeled 
by calculating the extreme velocities in each fiber and com- 
paring those to the expected velocities, is found to be lower 
than 3 km s" 1 at all radii. Subtracting the dominant sources of 
observational uncertainties from the observed dispersion in the 
rings, we find the intrinsic fiber-to-fiber line-of-sight dispersion 
is around 8 km s" 1 . 

Because the dispersion a is comparable to the rotation veloc- 
ity Vff, in the central regions, it may be necessary to correct the 
RC to obtain the circular velocity. If one assumes the disper- 
sion represents pressure support due to random motions, then, 
assuming axisymmetry and absence of radial motions, the cor- 
rection can be calculated and is found to be negligible (see e.g., 
S99). Alternatively, if one assumes the central regions are in 
a steady state, one can calculate the circular velocity v c using 
the virial theorem, 2T + H = W. Although the kinetic energy in 
rotation T and in random motions n depend on the distribution 
of the gas, and the gravitational potential W depends on the to- 
tal mass distribution, comparison to more detailed calculations 
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FIG. 2. — a) Hybrid rotation curve from SparsePak data (dots) and the 
HI data from S99 (crosses). The rotation curve from S99 is given by 
the solid line, b) Comparison of dBB's data (black dots) to SparsePak 
data within 2.5" of the major axis (grey dots), c) Comparison of dBB's 
and our rotation curve, coding as in panel b. 
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FIG. 3. — Radial velocity with respect to the systemic velocity, nor- 
malized by the rotation velocity, as a function of angle in the plane 
with respect to the minor axis. The dots, circles and triangles repre- 
sent the points at radii of 0.23, 0.45, and 0.77 kpc in the top panel, and 
1.71, 2.03, and 2.35 kpc in the top panel, respectively. The solid line 
represents the normalized curve for circular rotation. 



show the circular velocity is approximately equal to: 

v 2 c = v 2 v + f3a 2 . (1) 

Assuming an isotropic dispersion (/? = 3) results in a steeper RC 
(see Fig. 5). 

Because the origins of the dispersion and noncircular mo- 
tions are not known, and because the assumptions for the cor- 
rections may not be valid, these corrections are uncertain and 
merely give an indication of the possible range. 

3.3. Mass models 

The contribution of the stellar disk to the RC was calculated 
from the /?-band light profile presented in Swaters & Balcells 
(2002), and the contribution of HI was calculated from the HI 
surface density profile presented in Swaters et al. (2002), scaled 
up by a factor of 1.32 to account for helium. The stellar disk 
was assumed to have a vertical sech-squared distribution with a 
scale height zo = h/6 = 0.58 kpc, and the HI disk was assumed 
to be infinitely thin. For the dark matter halo we considered a 
generalized NFW halo of the form: 



HVJ ~ (r/r s r(\ + r/r s r^ 

where r s is the scale radius and po the central density. These two 
parameters are linked to the more commonly used parameters c 
and V200, depending on cosmology (see Section 4). 

We have fit mass models over a range of a and mass-to-light 
ratios (M/Ls). In each fit a and the M/Ls were kept fixed and 
the halo parameters were allowed to vary. Confidence levels 
were not calculated because the velocities and their errors are 
not free of systematic effects, and hence the \ 2 values were only 
used for relative comparison between the models. In Fig. 4 we 
plot the resulting \ 2 values as a function of a for an M/L of 1 . 
Fits to the uncorrected and pressure corrected RCs give virtu- 
ally identical results, with best fits in the range < a < 0.8. For 



the RC corrected using the virial theorem, best fits are found 
for somewhat steeper inner slopes, in the range 0.3 <, a < 1.0. 
For higher values of a the quality of the fits decreases rapidly. 

Best fitting mass models for an M/L of 1 are shown in Fig. 5, 
for a = 1 (the NFW profile) and a = 0. The difference between 
the a = and a = 1 models are predominantly in the inner 1 kpc, 
where the RC is most affected by noncircular and random mo- 
tions (Fig. 5). The best fitting parameters for a range of M/Ls 
and a = 1 are given in Table 1, both for the pressure corrected 
RC and the one corrected with Eq. 1 and (3 = 3. 

4. DISCUSSION AND CONCLUSIONS 

The high-resolution velocity field of the LSB galaxy DDO 39 
presented here reveals the presence of noncircular motions in 
the central regions and an intrinsic fiber-to-fiber dispersion of 
about 8 km s" 1 . The origin of either component is unclear; 
DDO 39 does not have a strong bar or spiral arms, nor does 
it appear to be interacting. Likely contributors are: star forma- 
tion activity, turbulence, and small scale structure in the disk 
or dark halo. Irrespective of their origin, the noncircular mo- 
tions and the large a/v^in the galaxy center make an accurate 
measurement of the circular velocity curve difficult. 

As a result, inner slopes in the range < a <, 1 are con- 
sistent with the data, and DDO 39 is compatible with a 
wide range of dark matter properties. Specifically, it ap- 
pears that the data presented here are consistent with the in- 
ner slopes expected in the currently popular ACDM paradigm. 
Although earlier simulations predicted inner slopes ranging 
from a = 1.5 (e.g., Fukushige & Makino 1997, Moore et 
al. 1999) to a = 1 (e.g., NFW), a recent study by Power 
et al. (2003) suggests that the differences in a are mainly 
due to resolution issues, and places an upper limit to a of 
1.2. Taylor & Navarro (2001), based on analytical argu- 
ments, find a may be as low as 0.75 at very small radii. 
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FIG. 4. — versus a, as determined from fits with M/L fixed to 
1, and the halo parameters left free. The full line gives the results 
for fits to the uncorrected rotation curve, the dotted line for pressure 
correction, and the dashed line for the virial correction of Eq. 1 . 
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Even though the data presented here eliminate several uncer- 
tainties that gave poor constraints on a from HI and Ha long- 
slit observations (van den Bosch et al. 2000, van den Bosch & 
Swaters 2001, SMvdBB, but see de Blok et al. 2001a,b), the in- 
trinsic noncircular and random motions still result in poor con- 
straints on a. Given that deviations from circular motions seem 
common in LSB galaxies (e.g., Walter & Brinks 1999, de Blok 
& Walter 2000), and are present even in a very LSB galaxy 
such as DDO 39, it seems that an accurate measurement of a 
and, from that, a tight limit on the nature of dark matter may 
be difficult to achieve observationally from a component that is 
susceptible to kinematic perturbations. 

Fortunately, the halo concentration parameter c is better con- 
strained observationally as it is derived from a fit to the RC as 
a whole rather than mostly from the central parts, although it 
may still be affected by systematic effects (SMvdBB). For the 
case of a = 1, the inferred halo parameters c and V200 can be 
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FIG. 5. — Mass models for the best fit NFW halo (solid line), for an 
NFW halo with c = 10 (short dashed line), and for a halo with a = 
(long dashed line). The dotted line represents the contribution of the 
HI, the dot dashed line the contribution of the stars with an M/L of 1. 
The contribution of the halo is similar to the best fit model, but is not 
shown. The black dots give the uncorrected rotation curve, the gray 
dots (only shown for the inner five points and offset by 0.03 for clarity) 
represent the rotation curve corrected using Eq. 1 assuming (3 = 3. 

compared to the values expected in a ACDM cosmology. Ex- 
pected values have a 2a range from 5 to 25, with an average of 
around 10 to 15 (NFW, Bullock et al. 2001). DDO 39's c ~ 5 
is at the low end of this distribution, and for different values of 
c the quality of the fit changes rapidly. For example, c = 10 is 
inconsistent with the data (see Fig. 5 and Table 1). 

The low value for c does not necessarily indicate an inconsis- 
tency with ACDM. For example, a bias toward low c values in 
LSB galaxies could be explained if LSB galaxies preferentially 
form in low density halos. In addition, the concentration param- 
eter also depends on the slope of the power spectrum of density 
fluctuations, and c ~ 5 is in agreement with models in which 
structure formation on small scales is suppressed (Zentner & 
Bullock 2002). Furthermore, c would be higher if DDO 39 
were closer than the adopted 12.8 Mpc. 

Accurate RCs determined from high resolution, two- 
dimensional velocity fields for a sample of LSB galaxies with 
intermediate inclinations, accurate centers, and regular kine- 
matics and morphology may provide RCs to measure c, or any 
parameter that depends on the global shape of the RC, with suf- 
ficient accuracy to provide useful constraints on the nature of 
dark matter. 

This research was supported in part from NSF grant AST- 
9970780. 
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